When metal organic frameworks turn into linear magnets 
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We investigate the existence of linear magnetism in the metal organic framework materials MOF- 
74-Fe, MOF-74-Co, and MOF-74-Ni, using first-principles density functional theory. MOF-74 dis- 
plays regular quasi-linear chains of open-shell transition metal atoms, which are well separated. Our 
results show that within these chains — for all three materials — ferromagnetic coupling of significant 
strength occurs. In addition, the coupling in-between chains is at least one order of magnitude 
smaller, making these materials almost perfect ID magnets at low temperature. The inter-chain 
coupling is found to be anti-ferromagnetic, in agreement with experiments. While some quasi-lD 
materials exist that exhibit linear magnetism — mostly complex oxides, polymers, and a few other 
rare materials — they are typically very difficult to synthesize. The significance of our finding is 
that MOF-74 is very easy to synthesize and it is likely the simplest realization of the ID Ising 
model in nature. MOF-74 could thus be used in future experiments to study ID magnetism at low 
temperature. 



PACS numbers: 75.10.Pq, 75.25.-j, 75.75.-c, 75.40.Cx 

The continued quest for the development of non- 
volatile memories and spintronic devices of smaller sizes 
requires the full comprehension of finite-size effects. To 
this end, over the last decade, exotic magnetic proper- 
ties have received much attention in experimental and 
theoretical studies-IS^SI Considerable emphasis has been 
given to the synthesis and prediction of materials show- 
ing mono-dimensional magnetism,'^!^^] qI^q referred to as 
ID or linear magnetism. While ID magnetism can be 
explained with the well-understood Ising model (dating 
back to 1925),'^ a satisfactory physical realization of 
this model in simple materials has not yet been found 
and I D magn etism is only observed in a few — often 
difScultPiSEll or dangerou^ to synthesize — synthetic in- 
organic materials and polymers. Although, for example, 
CrSb2 is one of the few materials that shows naturally 
ID anti-ferromagnetism, this property remains difficult 
to control and tune.'^^ In fact, theory has shown that 
strong spin fluctuations induce ferromagnetic disorder 
of ID-spin arrays at any temperature, independent of 
the extent of exchange interactions between neighboring 
gpj]-^g 115116] Thus, progress in the field of ID magnetism 
crucially depends on the availability of currently missing 
simple-to-synthesize model systems and materials. 

The main difficulty in engin eerin g good model systems 
exhibiting ID magnetism arePE^ (i) to find materials 
that have quasi- ID chains of spins with significant in- 
teractions and large magnetic anisotropy, (ii) to find ma- 
terials with a large ratio between intra and inter-chain 
magnetic interactions, (iii) to find materials in which fer- 
romagnetism is preserved at "reasonable" low temper- 
atures, (iv) to find materials with very few impurities, 
which tend to destroy ferromagnetism, and finally (v) to 
find materials that are simple, safe, inexpensive to syn- 
thesize, and where linear magnetism is easy to control. 
Historically, the engineering of ID magnetic materials 



has followed several routes. Attempts were made using 
inorganic materials such as Sr2Cu(P04)2, Sr2Cu03,'^and 
BaCo2V208,'^ along with non-periodic magnetic clusters 
or molecular magnets.'^ Another strategy is the combi- 
nation of organic molecules and tran sition m etals (TM) 
to form regular polymer ID magnets. f^ISEnilll xhe latter 
strategy offers a larger degree of freedom due to the high 
tunability of the diamagnetic organic separators , whic h 
can promote spin localization on the central TM.EEMl 

We propose here that metal organic frameworks 
(MOFs), a novel class of nano-porus materials, offer a 
versatile platform for the realization of ID magnets due 
to their high tailorability and tunability that results from 
their discrete molecular building-block nature. For 
this reason MOFs are already targeted in a large variety 
of applications such as gas-separation, gas-sensing, gas- 
capture, catalysis, and drug-delivery! ^^"^^ 1 In particular, 
in the following, we argue that the structural simplicity, 
low cost, and ease of synthesis of MOF-74 — together with 
the already existing understanding of this material — 
fulfill the criteria mentioned above and thus make it 
an outstanding candidate for studying linear magnetism. 
Note that signatures of ID-ferromagnetism in MOF-74- 
Co were already observed experimentally by Dietzel et 
al. in their pioneering work on this MOF.-^ From Fig. [l] 
it is apparent that MOF-74-TM (with TM = Mn, Fe, Co, 
Ni, and Cu) can be seen as ID magnet since it displays 
regular pseudo-chains of transition metals aligned along 
the basal plane. The helicoidal chains resulting from the 
atomic-motif of Fig.[T^ are interspaced by "long" organic 
linkers, suggesting that the inter-chain interactions are 
quenched. In fact, MOF-74 shows a large structural ratio 
(^ 3) between the separation of spins in a chain compared 
to chain separation (see Fig. [T]), establishing a required 
condition for the construction of a ID magnet. 

To elucidate the ID-like magnetic properties exhib- 
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TABLE I. MOF-74-TM net atomic charges (in units of the 
electronic charge), Qo and Qtm, and electron population of 
p and d orbitals on O and TM atoms, qo{2p) and qTM{3d), 
respectively. Magnetic moments, /i, are reported in units of 
Ms- 



TM 


Qo 


9o(2p) 


Qtm 


gTM(3d) 


M 


Fe 


-0.30 


4.75 


+0.50 


6.35 


3.625 


Co 


-0.95 


5.35 


+2.49 


5.17 


3.255 


Ni 


-0.61 


4.94 


+ 1.24 


8.33 


1.567 






FIG. 1. a) Frontal view of MOF-74, hehcoidal magnet 
chains are highlighted in green, b) side view of MOF-74, TM 
atoms are represented by green balls. dNN and dNNN are the 
nearest-neighbor and next-nearest-neighbor intra-chain dis- 
tances, while di-i is the inter-chain distance. The couplings 
^NN, ^NNN, and Ji_i, are defined in parallel. 



ited by MOF-74, we study the three isostructural 
materials MOF-74-Fe,l^ MOF-74-Co,l2^ and MOF-74- 
NiP To this end, we use density functional theory 
(DFT) with the PBE functional, as implemented in 
QUANTUMEsPRESSO.'2fil We employ ultrasoft pseudopo- 
tentials with wave-function and density cutoffs of 680 eV 
and 6800 eV. The pseudopotentials used for the TM (i.e. 
Fe, Co, and Ni) are also suitable for spin-orbit calcula- 
tions including relativistic corrections. The total energy 
is sampled with a 2x2x2 fc-point mesh, resulting in en- 
ergy differences converged to within less than 1 meV. 
Projected density of states onto selected atomic orbitals 
are performed on a denser fc-point mesh, i.e. 4x4x4. The 
SCF total energy convergence criterium is 1.4x10"^" eV. 
We need such tight criteria to be able to accurately sam- 
ple the delicate energy landscape originating from differ- 
ent spin arrangements. 

All calculations are performed on the experimental 
structures of MOF-74-Fe,'23lMOF-74-Co,'22!and MOF-74- 
Nijlll which crystallize in a rhombohedral primitive cell 
with 54 atoms and space group i?3. The calculation of 



the intra-chain J-coupling constants requires the freedom 
to have varying spin directions along a chain. But, the 
primitive cell of MOF-74-TM contains only 6 TM atoms 
that all belong to different chains (1 per chain), which 
does not give the required freedom. Thus, we construct 
a supercell extending the unit cell along the chain direc- 
tion, such that each unit cell now contains two chains 
with 6 TM atoms per chain, and a total of 108 atoms. 
Coordinates and relative lattice constants of the super- 
cells are reported in the Supplementary Information. 

Linear magnetism relies on ferro- or antiferro- 
magnetism that can only exist if the TM atoms have a 
non- negligible magnetic moment. We therefore begin by 
analyzing the localization of the magnetic moment on the 
TM atoms, combining the projected density of states and 
the Lodwin population analysis. The Lodwin analysis, 
similarly to the Mulliken analysis is an intuitive (but not 
unique) way of re-partitioning the electron charge density 
on each atom (and orbital), by projecting it onto indi- 
vidual orthonormalized atomic orbitals.^^ Table H] shows 
the Lodwin charges, relative contribution, and magnetic 
moments of the TM and O atoms in the three MOF-74 
investigated. The magnetic moments, fj,, were computed 
by integrating the spin-densities difference (pup — Pdown) 
of the d-p orbitals in the valence region of each TM. Al- 
though it is inadequate to draw decisive conclusions from 
the charge analysis of Table |Tj we observe that oxygen 
atoms in MOF-74 assume an interesting covalent nature, 
having repercussions on the final charges and magnetic 
moment of the TM in MOFs. We further confirm the 
local charge of Co in MOF-74-Co (+2.49), which was ex- 
perimentally assigned as 2+.!^ It is also interesting to see 
that the local charge of Fe in MOF-74-Fe behaves almost 
like the metallic case, thus increasing the local magnetic 
moment. The experimental magnetic moment for Co is 
4.67 Mb, which is larger than our computed value; a 
discrepancy connected to the well-known unphysical de- 
localization of the electron charge density that is intro- 
duced by the exchange-correlation functional adopted in 
DFT simulations Note that orbital magnetisnPM] is 
not included in our calculations, as its effect is typically 
very small.'^ 

Figure |2] shows the density of states of the valence 
bands projected onto the rf-orbitals of the TM atoms 
and p-orbitals of oxygen atoms (pDOS). Here we see 
that some of the electronic charge density of the TM 
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FIG. 2. Projected density of states onto Fe, Co, and Ni 
d-orbitals (gray) and O p-orbitals (red) of the valence bands 
of MOF-74-TM. Energy is given in eV with respect to the 
top of the valence band. Spin-up and spin-down densities are 
plotted above and below the zero line of each plot. 

spills-over (due to orbital hybridization) to the nearest- 
neighbor oxygen atoms. This diminishes the local mag- 
netization on spin-carriers and thus their total magnetic 
moment. Not surprisingly, the analysis of the pDOS to- 
gether with the charge analysis suggests that the mag- 
netization originates from the d-electrons (spin-down, 
see Fig. [2]) of the TM atoms. Note that the angle 
Ztm-O-tm ~ 90° ± 5° does not allow sufficient over- 
lap between the relevant orbitals enforcing the intra- 
chain ferromagnetism according to the Goodenough and 
Kanamori rules. '^'^ The above analysis clearly shows how 
the tunabihty of the organic linkers in MOFs can be uti- 
lized to increase the spin localization on the TM, and a 
more involved explanation can be found in Ref. 8, From 
this analysis we conclude the existence of localized mag- 
netic moments that can give rise to ferromagnetic cou- 
pling among TM atom chains. 

Although we have clarified the existence of chains of 
spin carriers, we still need to understand if spin chains are 
independent of each other (see Fig.[T]) in order to produce 
isolated spin arrays acting as linear magnets. To this end, 
a qualitative estimation of the magnetic independence of 
spin chains is obtained by performing calculations where 
each chain magnetization is assigned a random spatial 
starting orientation, which thereafter is free to relax to- 
wards the most favorably energetic orientation. If there 
is some degree of inter-chain interaction, each chain spin- 
magnetization will assume some preferred orientation. 
But, our results show that only a small rearrangement 
of the spin directions occurs, i.e. only ± 2° from the ini- 
tial directions, supporting the idea that chains are only 
weakly coupled. However, a quantitative measurement 
of such chain-chain interactions can only be obtained by 



TABLE II. Intra-chain J-coupling constants Jnn and Jnnn 
and inter-chain Ji_i for MOF-74-TM in cm~^. For clarity 
we report again the magnetic moment, n, in ^b, from Ta- 
ble|l] The standard deviation of Jnn and Ji_i is not reported 
because below the accuracy limit. 



TM 




Jnn 


Jnnn 


Ji-i 


Fe 


3.625 


28.1 ± 0.4 


6.0 


-1.2 


Co 


3.255 


40.1 ± 2.9 


4.9 


-1.9 


Ni 


1.567 


21.0 ± 3.5 


6.9 


-1.3 



calculating the inter-chain J-coupling constants, which 
follows next. 

Having established that the TM spin-carriers exhibit 
a substantial magnetization that can produce potential 
ferromagnetic coupling, our investigation moves to the 
calculation of the J-coupling interactions. Figure [T]d 
shows the magnetic pathways and defines the following J- 
couplings: The intra-chain Jnn and Jnnn, origin of the 
ID linear magnet properties; and, the unwanted inter- 
chain Ji_i interactions. A complete structural analy- 
sis shows that the intra-chain TM-TM distance falls be- 
tween 2.8 A and 3.0 A for MOF-74-TM, whereas the 
intra-chain distance falls between 7.5 A and 8.8 A, giv- 
ing reason to believe that the inter-chain J-coupling in- 
teractions are quenched. If each spin magnetization is 
constrained along the z-direction,'23 the coupling interac- 
tion, Jij , described by the complex Heisenberg-Dirac-van 
Vleck Hamiltonian simplifies to the ID Ising Modei^^ 

n 

H = -2J2J^JS!■S^ , (1) 

where is the projection of the spin operator along the 
z direction at site i. Due to the gyromagnetic factor, for 
the expectation values of S'f we use 1/2 of the magnetic 
moments /i in Table |Tj i.e. 0.813 for Fe, 1.628 for Co, and 
0.784 for Ni. We now use DFT to map the real system 
onto this model Hamiltonian by computing the energy 
differences of various ferro- and anti-ferromagnetic spin 
configurations, which in turn yields the J-couplings. Our 
supercell contains 6 independent TM atoms per chain 
(see Fig. [i]d) , resulting in 2^ = 64 possible different spin 
configurations, out of which only 16 combinations are 
linearly independent and compatible with our periodic 
boundary conditions. The coupling constants Jij are 
then obtained by solving an overdetermined system of 
16 equations with a least-square fit. Table [TT| reports our 
calculated values for the nearest-neighbor coupling Jnn, 
the next-nearest-neighbor coupling Jnnn, and the inter- 
chain coupHng Ji_i for MOF-74-Fe, MOF-74-Co, and 
MOF-74-Ni. Note that these calculations are a partic- 
ularly challenging task requiring high accuracy, as these 
energy differences are tiny compared to the total energy 
of a 108 atom unit cell. 

From Table [H] we see that the intra-chain J-couplings 
are larger and more positive than the inter-chain ones. 
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FIG. 3. Computed magnetic susceptibility xm (i^ cm^ 
mol~^) as a function of temperature T (in K) for MOF-74-Fe, 
MOF-74-Co, and MOF-74-Ni. The inset shows an enlarge- 
ment of the transition zone. 



suggesting the existence of linear ferromagnetism. On 
the other hand, the interaction among chains is very 
small and of anti-ferromagnetic nature. As expected, 
longer range J-coupling interactions, such as Jnnn, are 
of smaller magnitudes than the nearest-neighbor interac- 
tions and are expected to vanish at increasing distances. 
Although couplings for longer distances are in princi- 
ple easily obtainable from Eq. ([T]), such results are not 
presented here since they fall below our accuracy limit. 
Overall, the trend of the magnetic constants is main- 
tained between the three TMs. From our simulations 
the computed Jnn for MOF-74-Fe seems largely overes- 
timated from the experimental value of 4.12 cm~^, which 
was extrapolated by fitting experimental magnetic sus- 
ceptibility profiles.'SS On the other hand, our calculated 
inter-chain constant for MOF-74-Fe is in excellent agree- 
ment with the experimental result of —1.12 cm"^.-^ In 
summary, we conclude that the ferromagnetic intra-chain 
interactions are one order of magnitude larger than the 
anti-ferromagnetic inter-chain ones — confirming the pos- 
sibility of the existence of ID-magnetic phenomena at 
low temperature. 

We finally move to calculating the temperature- 
dependent magnetic susceptibility xm- Starting from our 
computed J-coupling constants, we can predict the mag- 
netic susceptibility, xm, which is measurable experimen- 
tally. We use Fisher's modeP^ 



XM 



Ng\ 



w(Jnn) = coth 



1 + m(Jnn) 
12/cfcT 1 — u( Jnn) ' 
/ kbT \ _ / 2hT 

\2Jnn/ V ^NN 



(2) 
(3) 



where N is the number of atoms in the chain, giso the 
g-factor, fcfc the Boltzman constant, and T the temper- 
ature. Figure [3] shows our calculated xm as a func- 
tion of temperature for the three MOF-74-TM investi- 
gated, using the JNN-coupling constants from Table In] 



The transition temperature corresponding to the phase 
transition from ferromagnetic order to anti-ferromagnetic 
order along the chains is given by the peak position of 
Xm- Obviously, the transition temperature depends on 
the J-coupling strength: the larger the J-coupling con- 
stant is, the broader the peak becomes and the higher the 
transition temperature. A similar dependence is found 
for the Xm magnitude itself, which decreases for increas- 
ing J-coupling constant. For MOF-74-Co, the tempera- 
ture dependence of xm was measured experimentally,'^ 
finding a transition temperature of 8-10 K, in good 
agreement with our calculated transition temperature 
of 13 K. The experimental maximum of the peak is at 
^0.17 cm'^ mol~^, whereas our calculated maximum is 
only at ~0.06 cm^ mol^^. However, this discrepancy is 
a result of the fact that our DFT calculated Co dipole 
moment of 3.255 fiB is too small compared to the experi- 
mental one of 4.67 /is (as mentioned above). ^"^ As can be 
seen from Eq. ^ the dipole moment /i enters the suscep- 
tibility as /i^. If we simply use the experimental dipole 
value, our peak maximum would be at ~0.13 cm"^ mol~^, 
in reasonable agreement with experiment. Furthermore, 
note that according to the susceptibility model used, the 
Jnn coupling constant for MOF-74-Fe has to be larger 
than the 4.12(6) cm~^ found experimentally through fit- 
ting data by Bloch et a/. such a small value results in a 
transition temperature too close to K and below the ex- 
perimental conditions reported in their study (2-300 K) . 
Equation ^ includes only the effect of Jnn, making this 
model quite unsatisfactory. The effect introduced by in- 
terchain Ji_i coupling constant in xm can be reintro- 
duced in Eq. ^ by replacing m( Jnn) with u(| Jnn/ Ji-i|), 
with the effect of slightly shifting all curves by ^ -3 K, 
bringing them in very good agreement with experimen- 
tal observation. Our estimated transition temperatures 
of all three investigated MOFs are clearly above the liq- 
uid He temperature, encouraging further experiments on 
linear magnetism phenomena in MOF-74. 

In summary, we have explored the existence of lin- 
ear magnetic phenomena in the metal organic frame- 
work materials MOF-74-Fe, MOF-74-Co, and MOF-74- 
Ni by using DFT calculations. Our results provide an 
understanding of the origins and magnitude of linear 
magnetic effects in these materials. We verify the exis- 
tence of intra-chain ferromagnetism and quenched anti- 
ferromagnetic coupling between chains, large enough to 
be observed at liquid He temperatures. The significance 
of our finding is that MOF-74 is easily synthesized, safe, 
and inexpensive. As such, it is likely to be the simplest 
realization of the ID Ising model in nature and has the 
potential to provide simple means to study linear mag- 
netism. Note that the MOFs studied here have not been 
tailored in any way to make them good ID magnets. In 
view of the high tailorability and tunability of MOFs, 
exciting new opportunities open up, where especially de- 
signed linkers and spin-centers decrease the spin-density 
delocalization and maximize the spin moments, resulting 
in larger J-couplings and higher transition temperatures. 
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